


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1975-06 


Control of nuclear propulsion plant power transients 


Chatham, Howard Maitland 


http://ndl.handle.net/10945/20793 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
3 Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
KN OX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


CONTROL OF NUCLEAR PROPULSION PLANT 
POWER TRANSIENTS 


Howard Maitland Chatham 








FCT TPS: 





T 
NV 


= 


POWER TRANSIENTS 
by 


howard Maitland Cnatham 


2 
CONTROL OF NUCLEAR PROPULSTI 
| 
| June 1975 


f 
f Meesis Acvisor: 





Approved for public release; distribution unlimited. 


T167532 








SECURITY CLASSIFICATION CF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE BERGE ECO Ee, ae 


REPORT NUMBER = F GOVT ACCESSION NO} 3. RECIPIENT'S CATALOC NUMBER 


ee 


4 TITLE (and Subiitie) 











5. TYPE OF REFORT & PERIOD COVERED 


Master's Thesis; June 1975 





Control of Nuclear Propulsion Pléent 


Power Transients 6. PERFORMING ORG. REPORT NUMBER 


7. AUTHOR(a#) 6. CONTRACT OR GRANT NUMBER(s3) 


Howard Maitland Chatham 


9. PERFORMING CRGANIZATION NAME ANO ADORESS 10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS | 


Naval Postgraduate School 
Monterey, California 93940 


11. CONTROLLING OFFICE NAME AND AOCDRESS 12. REPORT DATE 


Naval Postgraduate School June 197 
Monterey, California 93940 13. NUMBER OF PAGES 


4. MONITORING AGENCY NAME & ACORESS(i!f different from Controlling Oftice) 1S. SECURITY CLASS. fof thie report) 







Naval Postgraduate School Unclassified 
Monterey, California 93940 Se. Deere elcaoy DOWNGRAGING 






416. OISTRIBUTION STATEMENT (of thie Report) 


Approved for puolic release; distribution unlimited, 





17. DISTRIBUTION STATENENT (of the abstract snfezad in Glock 20, If different from Report) 


18. SUPPLEMENTARY NOTES 





18. KEY WORDS (Continue on twverse ice if necaacrer; and tdartilty by tock number) 


Nuclear propulsion, simulation, N.S. Savannah, power transients. 





20. ABSTRACT (Continue on reverse cide if necessary end identify dy block nimber) 


Tne nuclear propulsion pliant of the N.S. Savannah is simplified and 
Simulatec using Digital Simulation Language. The plant is subjected to an 
up power transient of sixty-five percent and a down power transient of eighty 
percent, vonstraints are imposed to limit power excursions, minimize reactor 
coolant hot leg temperature and. maximize boiler saturation temperature while 
maintaining avera:e coolant tempernzture eonstant, Combinations of control ro? 
movement and changes in reactor coolant flow rate are studied in order to 





ll A I AOS IS PETES OD, OIE ST sl 


DD yan 1473 EDITION OF 1 NOV 68 !E OBSOLETE 


S/N 0102-+014- 6601 —— eo eee Ta 
(Page 1) : SECURITY CLASSIFICATION OF THIS PAGS (When Data Entered) 





$2 CURITY CLASSIFICATION OF THIS PAGEMHher Dera Entared: 











+7 


determine the control systems that best satisfy the constraints and minimizs 
the time to steady state for each transient, 


NS SN cis CLS * OT LR NE TT 


Nha eens —wnerness Des: soe emp “Pes cE Spe: abe OE EO AT SEE TNS SG SRN: AAD A ASS a ee ES SG A i a ERE ek 


DD Form, 1473 





a Sem aol 


an fe 
S/N 0102-014-6601 SECURITY CLASSIFICATION OF THIS PAGE(When Data Entrro4) 








Control of Nuclear Propulsion Plant 
Power Transients 


by 


Howard Maitland Chatham 
Lieutenant, United States Navy 
B.S., Lowa State University, 1968 


Submitted in partial fulfillment of the 





“% 


by 
% 
+ 
Ole 


Thes 





thrar / 
V { 


ie) 


ABSTRACT 


The nuclear propuision plant of the N.S. Savannah is simplified 
and simulated using Digital Simulation Language. The plant is 
subjected to an up power transient of sixty-five percent and a 
down power transient of eighty percent. Constraints are imposed 
to limit power excursions, minimize reactor coolant hot leg tem- 
perature and maximize boiler saturation temperature while main- 
taining average coolant temperature constant. Combinations of 
control rod movement ana changes in reactor coolant flow rate are 
studied in order to determine the control systems that best satisfy 
the constraints and minimize the time to steady state for each 


transient. 
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I. INTRODUCTION 


As the energy crisis becomes more pressing every day, nuclear 
power stands out as a most promising solution. Generation of electric 
power from nuclear energy was first accomplished in 1951 at the Atomic 
Inery Commission's (AEC) test facility in Idaho. Since then nuclear 
power stations have been built for commercial use in many countries 
and it has been estimated that three quarters of the total generating 
capacity of the United States will be from nuclear power by 2000 AD. 

But power generation on land at present takes a second seat to 
marine plants. The United States alone has over one hundred ships 
with nuclear propulsion systems. mussia, England and France also 
have vessels with similar propulsion. The initial construction 
costs have limited most of these to naval craft, though West Germany 
and Japan have developed commercial ships with nuclear propulsion. 
The United States also built a&@ commercial nuclear vessel, the Nuclear 
Ship (N.S.) Savannah, which has since been decommissioned. 

The basic advantages of nuclear power versus fossil fuels are 
the small amounts of fuel required for substantial power generation, 
the relatively few pollutants generated, and the longevity of the 
fuel before refueling is necessary. Initiai set up costs and costs 
of refueling are considerably more than for fossil fuel plants 
but the fewer number of plants for tne same power generation and 
the longer lifetime of the plants more than counter-balance these 


costs./ 8/ 


Marine nuclear power plants differ from stationery power generating 
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plants. Normally they are smaller physically and the output size 
required is less. In addition shipboard plants have more substantial 
power transients imposed on them. an electric generating station 
will have periods of peak load changes from normal steady state 
demands but these are insignificant compared to the propulsion 

orders that can arise in a maneuvering situation at sea. Power 
transients of up to eighty percent are not unlikely and a different 
mode of control than used on stationery plants is required. 

Control systems on a nuclear plant keep system variables within 
proper limits and shut down the reactor when possible damage to the 
reactor or the environment is about to occur. Though the pressurized 
water reactor, the normal plant for shipboard use, is inherently 
stable it is also quick to respond. An automatic contrel system 
for this propulsion plant would minimize the possibility of human 
error, be able to cope with the quickness of response, and could 
also reduce the number of watchstanders required to operate the 
plant. 

This thesis investigates the automatic control of a nuclear 
propulsion system to large power transients. To insure the validity 
of the system, the plant is modelec as closely as possible to that 
of the N.S. Savannah, the only unclassified power plant for which 
sufficient data was available. Similar load changes as applied 
in this study were conducted during power range testing on the 
N.S. Savannah with satisfactory results./ 9/ 

The thesis is organized into sections as follows: first, the 
N.S. Savannah is described paying particular attention to the pro- 


pulsion systems; next the reactor kinetics and thermodynamics equations 
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for the plant are developed. Following these are sections covering 
the results of the power transients to the plant without control 
systems and then with various control schemes applied. Finally 
the conclusions ana recommendations for further study are discussed. 
An appendix follows presenting the comments and problems on the 


computer simulation. 
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II. N.S. SAVANNAH 


Pee Lois SHIP 

The Savannah was designed as a passenger-cargo ship and as 
such even with & conventional propulsion plant might not have been 
economical compared to bulk cargo ships. She was built to promote 
the acceptance of nuciear vessels throughout the world and to serve 
as a prototype for future merchant vessels. Because of the economical 
problems she was subsidized by the government but run by civilian 
concerns. 

The ship began operations in early 1962. She made several trips 
overseas but political pressures in several countries against 
nuclear power curtailed many port visits. After much "showing 
of the flag" at home and abroad she was decommissioned in late 1971. 
Even though economically unsound at the time the ship was termed 
a success due to the performance of her propulsion plant, proving 


that commercial nuclear systems were practical./ 1/ 


B. THE POWER PLANT 
Tne basic nuclear power plant is divided into two system, the 
primary and the seconaary. 


1 Tre Primary system 





This system generates the heat to boil water in the secondary 
system. It consists of the nuclear reactor, coolant pumps and 
piping, the heat exchanger and auxiliary systems. All of these 
are enclosed in a containment device which prevents tre release of 


radioactive particles in case of some nuclear accident. 


16 


| 
‘ 
*. id 
* « 


“ 


- . a 
! “ 
sf 
‘ . 
7 oe 
® 
° * ‘7. > 





The reactor core is made up cf uranium dioxide formed into 
thin fuel plates coated with a cladding material. ‘Twenty-one control 
rods made of enriched boron-stainless steel are inserted between 
groups of fuel elements. The rods absorb neutrons and when fully 
inserted into the reactor sufficient neutrons are absorbed to prevent 
a self-sustaining nuclear reaction. As the rods are withdrawn less 
neutrons are absorbed tnus allowing more fissions to take place. 

The fission process generates thermal energy which is con- 
ducted through the cladding to water circulating through the reactor 
core. The water is kept pressurized to prevent boiling in the 
system. Tne water is forced from the reactor by coolant pumps in 
two loops. Midway in each lcop is a counter-flow heat exchanger. 
This method of heat transfer to the secondary system is used to 
insure there is no possible contamination of the propulsion system 
by nuclear by-products in the coolant. 

The major auxiliary systems consist of a pressurizer, coolant 
purification and make up systems. A steam buble is created by 
electric heaters in a pressure vessel pipec to the coolant loop. 

The higher temperature in the pressurizer maintains the bubble in 

this system and not in the primary loop. The bubble allows for 
thermal expansion and contraction in the primary system while 
maintaining a relatively constant pressure. The purification system 
draws off twenty gallons perminute of coolant to run through an 

jon exchanger. The make up system stores the purified coolant and 
charges additional water as necessary to the primary loops to maintain 
volume ae), 

Data for analysis of this and the secondary system are 


included as Table 1./8,97 


le 








TABLE 1 


Reactor maximum operating power 
Reactor flow (total) 

r. (maximum operating power) 

T. (maximum operating power) 
eave 

T. (maximum operating power) 

ue (maximum operating power) 
High power scram 

High qT scram 

Steam flow (total) 


Rod speed (maximum) 


Fuel loading 
Ee 


y278 


Volume fractions 
water 
Control rods 
fuel 


stainless steel 


Ap 


18 


N.S. SAVANNAH PRIMARY AND SECONDARY CHARACTERISTICS 


70 Megawatts 
8.64x10° lb/hr 
520°F 

496 °F 

508° F 

1108°F 

45a. EF 

130% 

540°F 

2.67x10° lb/hr 


20 inches/minute 


312.4 kg 


Ciao ke 


Oe>bT 

0.041 

0.247 

0.145 
-2.3x107?/°F 
-1.88x1077/°P 
715x107" 

5x10"? seconds 


0.1/seconds 








2. The Secondary System 


This system uses the heat of the primary system to create 
steam in a boiler connected to the heat exchanger. The propulsion 
system is very similar to most steam turbine systems. Steam from 
each boiler is piped to a propulsion turbine, a turbo-generator, 
the steam-driven feed pumps and auxiliary systems. The propulsion 
turbines also have a steam dump which allows steam to be by-passed 
around the turbine to the condenser. All condensed steam is pumped 
back to the boilers by feed pumps. 

4. The Plant Model 

Due to the complexity of a nuclear power plant some simpli- 
fication of the model was undertaken in order to carry out this 
study. A single coolant locp primary plant was assumed. All primary 
auxiliary system effects were neglected. Complete mixing in the 
reactor and a constant differential rod worth were also assumed. 
The ship's electrical loads (hotel loads) and auxiliary steam loads 
were assumed to be constant regardless of the power level. Feed 
water temperature and flow rate were also assumed to be constant. 

The actual variable available to the watchstanders on the 
Savannah were steam flow (m_), coolant flow (m,), boiler pressure 
and temperature CaeeG coolant hot leg, cold leg and average 
me anata 3 temverature differential across the 


heat exchanger (4T), neutron power (n), and control rod height./ 8/ 


temperatures (T 


In a saturated steam system P. corresponds to Te and this variable 
was not used in the simulation. As the reactor thermal power is 
equal to the temperature differential across the heat exchanger 


multiplied by the coolant specific heat and the coolant flow rate, 
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me andAT were combined into one term which when normalized became 
P. Normalization cf steam flow and neutron and thermal power allowed 
comparison of the three power levels on the same plot. Other plots 
generated were control rod height and system temperatures. Control 
rod height was zero referenced at one hundred percent steady state 
power and represented the height of the group of rods which was 
controlling the reactor. 

Figure 2-1 shows a simplified diagram of the power plant 


model. 


- 
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ITI. REACTOR KINETICS AND THERMODYNAMICS 


A. THE REACTOR KINETICS EQUATIONS 

Reactor kinetics is the study of the time dependent behavior 
of a nuclear reactor. The first equation equates the rate of change 
of neutrons causing thermal fission to the rate of production of 
prompt, delayed and source neutrons causing thermal fission less 
the rate of loss of neutrons of all types. The other six equations 
relate the rate of change of the six delayed neutron precursors 
to the rate of production of each group of delayed neutron precursors 
less their rate of loss. 

To simplify the reactor model the six groups of delayed neutron 
precursors are lumved into a sinyle group./ 5/ The decay constant 


becomes dependent on the variation of Ke from one, defined as 


a 
reactivity, but has been assumed constant for this study with negli- 
gible effect on the results. In addition, since the contribution 


of source neutrons while operating above one percent power is insig- 


nificant, tnis term was also neglected. Thus the equations are as 


follows: 
n=n1-@8)+Ac-_n (3-1) 
« x 
IK ee 
¢ =n,.B-A-c (3-2) 
th 


For use in the computer simulation and power comparisons, 
normalization of these equations is convenient. oy These normal- 
ized forms are shown on the next page with N=n/n(0) and C=c/c(0). 


The initial quantity is referenced to one hundred percent power. 


ee 





Ni =a arate» S atnines Gn (Onan ; 
B(heeg = 1) + (C= ¥) (3-3) 
eit. 
C= A(N - C) (3-4) 
Kore is the ratio of the number of neutrons produced from fission 
in one generation to the preceeding generation. When K equals 


eff 
one the reactor 1s said to be critical and neutron population is 
constant for successive generations. The power level is also constant 
and a steady state condition exists. Increasing Kore causes increased 


neutron production, the supercritical condition. Decreasing Kore? 
Subcritical, causes a decrease in neutron production./ 7/ 


Three major factors affect Ke These are average coolant 


er 
temperature in the core, average fuel temperature and control rod 
height./11/ An increase in coolant temperature decreases the density 
of the water which causes less moderation, slowing down of fast 
neutrons, and fewer thermal fissions occur. Thus Kore 1s reduced. 
Numerically this is computed by multiplying the change in temperature 
by a negative constant called the temperature coefficient of reac- 
tiveicy, in this case X,. A similar effect takes place due to the 
temperature of the fuel and this negative reactivity coefficient 
is Ape 

The height of the control rods also changes reactivity. As 
rods are withdrawn more reactivity is "inserted" in the core; driving 
rods in reduces the reactivity. The amount of reactivity per inch 
of travel is called differential rod worth (DRW). For the Savannah 
the DRW was computed to be 7.0x107*/ineh./ 8/ The zero reference 
for rod height was set for one hundred percent steady state power 


and the control rod worth (CRW) is found by multiplying the Daw 


by the rod height. 


20 





Combining all these terms and noting that -0.120988 is the 
reactivity due to Tat 508°F and T, at 1108° F, the following 
equation results: 


Kopp = 10120988 +a ,-T toe: T, + CRW (Gas) 


B. THE THERMODYNAMIC EQUATIONS 
1. Reactor Heat Transfer Equations 
The total power generated in the reactor equals the thermal 
capacity of the fuel times the rate of change of the fuel temperature 
plus the heat transfer coefficient of the fuel times the temperature 
difference between the fuel and the coolant in the core. 


& 


4 = mee wee - h(t, ~ Dee) (3-6) 


And the latter term equals the thermal capacity of the coolant 
times the rate of change of the average coolant temperature in the 


core plus the coolant flow rate times the specific heat of the coolant 


times the temperature difference of the coolant across the reactor. 


€ 


= . a . - se 
Ne ae acon, ins “rp (Toh Tee? (3-7) 


+ 


hi (Tp ~ “Cav 


where T= ae + Bane (3-8 ) 
2. Heat Exchanger Itquations 

The coolant flow rate times the specific heat of the coolant 
times the temperature difference of the coolant across the heat 
exchanger equals the thermal capacity of the coolant times the 
rate of change of the average temperature of the coolant in the 
heat exchanger plus the heat transfer coefficient of the heat ex- 
changer times the temperature difference between the coolant average 
temperature in the heat exchanger and the saturated liquid in the 


bowler, 


iy) (3-9) 


oa chs (T, - es “op eae are co 
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And the latter term is equal to the thermal capacity of the boiler 
liquid times the rate of change of the saturation temperature in 
the boiler plus the steam flow rate times the difference of enthalpy 
between the steam and the feed water. 
a eens as Oe Ty : ce coins (3-10) 
eo = a ~ ye Co) 
3. Modification of the Point Fquations 
The above equations, representing the point nodel Ji 
do not account for the transport delays of six seconds in the coolant 
ieep nor in the boiler. Letting 
T ant?) =: T, (+6 ) (a2) 
and T(t) = TL, (+6) a5) 
the time rate of change of Ty and a can be evaluated by expansion 
by a Taylor series and neglecting second order and higher derivatives. 
In addition the enthalpy difference in equation (3-10) is approx- 
imately constant for this plant. 
Evaluating constants as necessary, including coolant flow 
effects, and normalizing all power and flow terns / 3/, the following 


equations were used for the simulation: 


oe (3-14) 


= (600 N - oe 


: Cav 

ena : (T - ene eet cae / TQ) 1/25 Cai 
ee el (ee oe eG (3-16) 
Toe ate Dy) © (3-17) 
= 70 (Pp -¥)/3/M (3-18) 

: -19) 

aee (Ty, i T,)/2 (3-19) 
TL. ((2 - 0.48/M) T, + 0.96 t /i)/(2 OAS hy: 5220) 
P = 0.02 (Unie -T,). e721) 
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IV. TRANSIENT ANALYSIS WITHOUT CONTROLS 


A. TRANSIENT SELECTION 

The transients selected for this study were an up power maneuver 
of twenty to eighty-five percent and a down power maneuver of one 
hundred to twenty percent. Both these transients were run on the 
N.S. Savannah and its control system handled both these transients, 
except as noted in section VI, satisfactorily./ 8/ The up power 
maneuver represents a load increse from normal "hotel loads" to a 
"full bell" (normal full speed) in ten seconds. The down power 
test represents the tripping of the propulsion throttle from maximum 
power, a three second transient. 

In the up power case steady state conditions were assumed to 
have been reached when Ne was within three degrees of its finai 
value and neutron power within four percent power of its final value. 
For the down power case the temperature requirement remained the 
Same but due to the more significant overshoots the power band was 
reduced to two percent. 

It should be noted that due to the relatively low fuel enrichment 


of this reactor a significant ~%, is generated. Without this term 


if 
Tie. would have to return to its initial steady state value after 
any transient if no controls were applied. This characteristic 

is used in many power plants where more highly refined uranium can 
be used. In the Savannah this was not feasible but a constant 

Eee control was still desired thus rod movement was used to counter- 


act the fuel temperature effects. /13/ Subsequent sections discuss 


the restrictions and advantages of the ar control system. 
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Bereta UP POWER TRANSIENT 
The up power transient was rin without any control system to 
cetermine the natural system response. Steady state was reached 
in ninety seconds after the start of the maneuver. As anticipated, 
as T, increases with power, to maintain K 


it eff 


decrease. AT increases as does the temperature difference between 


equal to one T must 
ave 


m m + zl i 
es and te! Figures 4-1 and 4-2 present the plots of Tae 


Th and a and Y, N and P versus time. 


Steady state values were as follows: 


TF foe ~ = 85.0% 
Tove = A68.5°F N = 81.3% 
= A5G.5 © P = 84.6% 
TO = APon ene 


The peak overshoot of N was nine percent power. The minimum 
ue was 424.0 F. It wes noted that Ts was the first temperature to 


respond, followed in order by To ane and T This is due to the 


ne 
sudden change in steam demand draws energy from the boiler which 
reduces Tse The water in the heat exchanger has more heat removed 
thus cm decrease. es necessarily decreases and since the thermal 
power of the reactor is less than the steam demand T, increases 


h 


relative to Toe But due to the rapid decrease in TO qT decreases 
also. Also noted was the fact that neutron power lagged the thermal 


power in the initial response by five seconds and it was over thirty 


seconds before it caught up with P. 


C. THE DOWN POWER TRANSIENT 
The natural response of this transient reached steady state in 


seventy seconds. Just the opposite of the up power transient occurred 


ah 
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in the direction of tne temperature values though the order of 
response stayed the same. As before the neutron power lagged thermal 
power. The initial delay was about three seconds and the crossover 
point was not reached until seventy seconds. The peak overshoot 
was only 2.3% power. (This represents a 11.575 overshoot of the 
steady state value as compared to only a 10.675 overshoot on the 
up power transient.) 

Figures 4-3 and 4-4 show the results of this transient versus 


time. Steady state values were: 


T, = 559.8°F ff = 20.0% 

eee N = 241.2% 

T. = 554.8°F P = 21.1% 
= oT 

i DAGne as 


The maximum Ty Wagoe. fH. 
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Ve. CONTROL FOR THE UP POWER TRANSIENT 


A. BASIC CONTROL 

The normal control system for a pressurized water reactor is 
to keep ee in a specified band. The primary reason for maintaining 
Caan Im this transient is to raise ea Without controls it was 
noted that un dropped to 423° FP, The higher Te can be kept the more 
efficient the secondary plant will be. The Savannah used a control 
of TY, = 508 4 3°F which was adopted in this study also. Since 
in this case Pew decreased rod withdrawal was required to add 
reactivity to the core and raise system temperatures. The maximum 


rod withdrawal rate of twenty inches per minute was used to help 


speed up the response of the plant. 


B. CONTROL LIMITS 

In order to prevent excessive power generation in the core from 
causing a reactor accident, an automatic insertion of ali control 
rods (scram) is initiated when N reaches 130%. Due to the tolerances 
involved in the power detection instruments and in the scram cir- 
cuitry a ten percent power band below the scram set point was 
determined to be the maximum power allowable in this simulation. 

To insure that this limitation would not be exceeded by any 
transient a test was run using a twenty to one hundred percent 
maneuver. When N reached 120% rod motion was terminated but power 
kept on increasing until it just reached 130%. It was evident 
that a lower power limit was required. Using a 1105 limit the max- 


imum power attained was 117,53. Limits in between the two limits 


a 


either exceeded or were extremely close to 120% thus the 110% was 
adopted for further transients. 

Attempting the test transient with this limit resulted ina 
maximum power of 118% but a further difficulty was encountered. 
Since rod motion can only counter less than 1.25°F per second, 
Tove had originally dropped but was increasing and had reached 
495°F when rod travel ceased. 8 continued to increase and when 
power dropped below 11075 rods were bumped to maintain that level. 
When ess reached its lower limit of 505°F rod motion again ceased 
but as lee rose above 511° F rod insertion commenced. By the time 
Tete turned at 513°F power had been driven down to 59% and was 
still decreasing. 

This over action necessitated a rod position limit. The rod 
worth corresponding to T|, at 508°F and dependent on steam demand 
was used as the limit. Since T. at steady state equals Tove plus 
Y times 600, the limit was 0.0138-:(yY - 1). This corresponds to a 
rod height of 19.7:(¥ - 1) inches. Not allowing rod withdrawal 


past this point precluded the extreme power and temperature 


oscillations encountered above. 


fee ons ROL SYSTEMS 
in. Control 
ave ————— 
The 20 to 85% transient was run using the above limits and 


steady state was reached in 150 seconds. Steady state values were: 


T = 515.5 F Y= 85.0% 
ee OO) N = 87.9% 
T. = 494.7 F P = 85.3% 
tad £62.35 2. 


D4 


. 


eae 





Figures 5-1, 2 and 3 show the plots of temperatures, power 
levels and control rod height versus time. The minimum ae reached 
was 443.2°F, well above the minimum for the system response without 
controls. Maximum power was 118.2% and though N again lagged P 
initially, the crossover point was reached in only eighteen seconds, 
substantially less than in the previous run. 

It was also noted that N was still in excess of ¥and P 
after 200 seconds. This was due to the fact that the fuel with its 
slower time response than the coolant had not reached its final 
value. In order to raise T, the neutron power must be greater than 


i 


the thermal power anc as the steady state value of 7, 1S approached 


f 
the difference in power levels will diminish. 

Initial rod motion was at the maximum rate of twenty inches 
per minute. But during the period of bumping rods to maintain the 
110% N limit an average rate of 8.67 inches per minute was obtained. 

Zee aniticivatory Control 

In an effort to reduce the response time, early rod with- 
drawal was initiated. As soon as steam demand exceeded neutron 
power by a five percent tolerance band the rod pull was commenced 
instead of waiting until Len exceededits band. 

In this case steady state was reached in 135 seconds even 


though rod travel was only initiated five seconds sooner than in 


the previous case. Steady state values were as follows: 


a 515.2 F Y = 85.0% 
© 4 
= e =- 87. ) 
ieee 505.0°F N Loe 
Pert 494.8°FR P = 85.4% 
— ¢ .5°F, 
Ue = 462.5°F 
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Minimum Ts was only 446°F in this case, an increase above 
the previous run. Maximum power reached 119% and the crossover 
point for N and P was only fifteen seconds. Again the initial 
rod withdrawal was at twenty inches per minute whereas the subsequent 
average rate was 9.12 inches per minute. See figures 5-4,5 and 6. 

3. Discrete Flow Change 

Many reactor plants have the capability of two speed coolant 
pumps. At high power levels, normally above 503, the pumps are 
shifted to the faster speed usually twice the flow rate of the 
slow speed. This reduces the AT and with ae constant decreases 


iT Meets 16 a desiriqus result as wil! be discussed {or tae 


h° 
down power transient, an investigation of this flow change was 
conducted for the up power maneuver too. 

A simulation was run without temperature control but with 
the coolant pump shift at 50% N and steady state was reached in 
ninety-four seconds, only slightly slower than without the pump 
shift. 

Another run using the a control and the pump shift reached 
steady state in 156 seconds only six seconds slower than without 
the flow change. 
Combining all three control schemes steady state was reached 


in 141 seconds. This also was only six seconds slower than the 


anticipatory control case. Steady state values in this case were: 


1 = 510.1 F ‘= 85.04 
f _ e a i 
a 505.0 F N = 87.9% 
ae 499.9°F P = 85.8% 
e =TAGee ee 
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The system showed an increased Te minimum of 447.1°F and 
a reduced power level maximum of 117%. More of an overshoot in P 
was observed than before, 2.5% power, and the second rod withdrawal 
average rate was reduced to 8.66 inches per minute. Figures 5-7, 
8 and 9 show the pertinent data. 

4. Continuous Flow Change 

in an attempt to reduce the response time a continuaus flow 
change was assumed. The coolant pump was assumed to be variable 
speed, normal flow rate at 0% power and twice that rate at 100% 
steam flow. In between a linear relationship was assumed thus 
M equaled ¥Y plus one. Steam flow was used in this control because 
N was known to exceed 100% whereas steam flow was limited by the 
boiler and turbine construction. 

In this case the transient system reached steady state in 
138 seconds, three seconds faster than in the previous case. 
Minimum T_ was 447.8 F and maximum power was 117’. Figures 5-10, 


11 and 12 shcw tre variable plots and steady state values were: 


tT = 510.5° F Y = 85.0% 
6 at 
= e F Y = 8 .o% 
es I 78; 
rT 499.5°F P = 85.7% 
mo _ } 
Tr. = 462.1°F. 


The P, N crossover point was just over fifteen seconds and 


the second rod withdrawal average rate wes 8.76 inches per minute. 
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VI. CONTROL FOR THE DOWN POWER TRANSIENT 


Pema LG CONTROL 

As with the up power transient the control for the down power 
transient is based on maintaining ae = 508 + 3°R, But instead 
of attempting to raise Tos decreasing Ty is one of the primary 


concerns. High T, could cause fuel cladding burnout and release 


h 
of fission products to the coolant. The N.S. Savannah had a high 
Th scram at 540 F and this limit was adopted for this study. As 
temperature indicators are more accurate than power instruments 

the large tolerance band needed for the high power scram is not 
needed for the high Ty scram. 

Another concern is power overshoot, 1.e., power decreasing 
below the twenty vercent steam demand. Many reactor accidents are 
more severe at low power levels and levels of less than one percent 
virtually shut down the reactor. 


Since oe needed to be decreased, as seen in figure 4-4, rods 


were inserted into the core. 


B. CONTROL PROBLEMS 

eee control was attempted using several rod speeds and control 
limits. Using the maximum speed power decrease to less than one 
percent. A minimum power limit of 15% was established to allow 
rod motion after N had past Y thus continuing to reduce Tove but 
keeping power above several percent. To preclude the temperature 


and power oScillations noted in the up power transients the same 


rod limitation was imposed. Thus rods could not be inserted past 
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19.7 (¥ - 1) inches, the steady state height for the rods when 
{ was 508°F. 
ave 


Even using these limits either the reactor reached power levels 


less than one percent or T, exceeded 540° F. This was also the 


h 
case on board the Savannah. To counteract this situation an automatic 
steam dump was designed. As steam pressure suddenly increased, as 

in a down power transient, the steam dump was automatically opened 
fully into the condenser. ‘This throttle then could be shut at a 
slower rate while allowing the propulsion turbine to maintain its 
proper load./ 8/ 

No mention of a specified cinsure rate of the steam dump throttle 
could be found in the literature. Various length transients were 
Simulated at the slowest rad speed, five inches per minute. A 
sixty second transient was selected as being easily handled mechan- 
ically though a maximum ty od 550°R was obtained, 

Using this maneuver several rod speeds were tried to decrease 
Ty without shutting down the reactor by too rapid insertion of 
reactivity. Ten inches per minutes was found to reduce Ty just 


below 540 PF while keeping power over seven percent. This rate 


was used in all down power controls. 


C. THE NEW TRANSIE 

The sixty second down power transient was run without control 
systems to establish reference data. Steady state was reached 
in 108 seconds. The variable plots are shown in figures 6-1 and 
6-2. As noted in the up power transients neutron power lagged 
thermal power initially. The crossover point was reached in 115 


seconds with an insignificant neutron power overshoot. The steady 
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State values of temperature remained the same as for the three 


second transient while N was 22.0% and P was 21.0%. 


Dee CONTROL SYSTEMS 
te Control 
ave ———————— 
With this control system the steady state conditicns were 


reached after 356 seconds. Figures 6-3, 4 and 5 show the pertinent 


plots. Steady state values were: 


T = 510.4°F Ty =) 20n 0% 
iP = 508.0°F N = 18.0% 
ave 

TOO 505.6°R P = 19.8% 
T = 498.1°F, 


Maximum Th was just under 540 F and minimum N was 7.4%. 

N reached P after about thirty seconds and maintained about the 
same rate of decrease until thermal power reached 22%. Because 
of the longer and less peaked power overshoot the subsequent rod 
motion was not linear at first. After twenty seconds an average 
rate of 1.63 inches per minute was obtained while the initial rod 
travel had been at a rate of ten inches per minute. 

It should also be noted that the neutron power exhibits oscil- 
lations of up to two percent power even after steady state. Though 


L i ; To and Te have stabilized T, 1s stiil decreasing and 


h’ “ave 
as N is affected the most by Tp and 1S more responsive to reactivity 
changes at low powers these oscillations will coniinue for several 
minutes. 
2. Anticipatory Control 

° 


In this case rod insertion was commenced wnen steam demand 


was less than neutron power by five percent. The response time 
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was greatly reduced. Steady state was reached in 315 seconds, over 
forty seconds faster than with the Toye control only. Figures 6-6, 


{ and 8 show the relevant plots. Steady state values were: 


TOF 510.5 F Y = 20.0% 

Tive = 208.1 F N = 18.0% 

TOF 505.7 F P = 19.8% 
_ 0 

eS Gera 


Th maximum was reduced to 536.5 F and minimum power was 
increased to 7.7/5. Rod insertion began with in four seconds. 

N decreased at a more rapid rate than P for most of the steam demand 
change. 

The rate change in the initial rad response is due to rod 
initiation at 95% steam demand but since N followed ¥ so closely, 
just five percent or less difference until steam demand stabilized, 
only a 6.05 inches per minute average rate was abtained until oa, 
reached 511 F. Even so, 9.31x107% p was inserted before rod motion 
would have been started in the een control mode. As with the ee 
control the subsequent rod motion was not linear until well after 
N returned to 1555 and then maintained a rate of 1.77 inches per 
minute. 

3. Discrete Flow Change 
By using a pump shift when increasing power the reactor 
T is reduced. Thus ee is reduced when operating at high power 
levels. Down shifting pumps when power decrease to 50% was attempted 
and the response time was further reduced to 309 seconds. Figures 


6-9, 10 and 11 show the power, temperature and control rod height 


oO = i ° 
plots. The maximum Ty was 535.8 F and power reached i. 7% again. 
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Steady state values were: 


Is = 510.4°R W = 20.0% 
T = 508.1°F N = 18.0% 
ave 

we = 50507 e P = 19.87% 
T, = 498.1°F. 


In this instance, N exceeded steam demand slightly more 
than in the previous case and the average initial control rod rate 
was 7.4 inches per minute until ie reached 511° F. Thus 15.6x1074 
f was inserted before the ten inches per minute rate was established. 
The second rod motion when in the linear area was at a rate of 
1.59 inches per minute. 

4. Continuous Flow Change 

As in the up power case a flow rate of M=1+¥ was 

established and the results are shown in figures 6-12, 13 and 14. 


Oteady state was reached in 310 seconds and the values were: 


= 510.0°F Y = 20.0% 
7 = 508.0°F N = 18.0% 
ave 

TS 506.1°F P = 19.7% 
T = 498.1°F. 


Ty reached a maximum of 535.3 F and power again had a 


minimum of 7.7:3. The control rod travel rates where insignificantly 


Gifferent from the discrete pump shift case. 
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VII. CONCLUSIONS AND RSCOMMENDATIONS 


poe CONCLUSIONS 
1. The Up Power Transient 
In this case several specifics were noted. The sooner the 
initial rod withdrawal, the quicker was the response of the- system. 
Increasing the coolant flow rate did not substantially degrade the 
response time but did reduce the maximum power and increased the 
mininun TS: The continuous flow change system gave a slightly 
better response than the discrete flow change. 
2. The Down Power Transient 
As noted above the sooner the rod motion started the quicker 
the response. But in this case a ciscrete flow change further 


decreased the response time as well as the maximum T The minimum 


h’ 
power level was not affected. The continuous flow change caused a 


slight decrease in T, while response time increased insignificantly. 


h 
foe banal Control 

A satisfactory control for these transients was developed 

using only two rod speeds, twenty and ten inches per minute. A 

simple two speed motor for the control rod drive mechanisms could 

be used, greatly simplifying the variable speed used on the Savannah. 
Overall comparison showed that the control using the — 

band, an anticipatory rod motion response and continuous flow changes 

was the system which responded best to both the up power and the 

down power transients. But due to the complexity of variable speed 


motors vice dual speed the discrete flow change contro! would be 


amore practical system with only a slight reduction in response. 
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Be. XECOMMENDATIONS FOR FURTHER STUDY 

Several extensions of this investigation merit further study. 

1. This study used a two speed rod control. Either a variable 
Speed system as used on the N.S. Savannen or optimizing the rod 
speed to one value, thus enabling the use of a still simpler motor, 
could be used as a2 basis for further work. 

2. As this study looked at the two transients separately, a 
study of the two maneuvers in rapid succession would provide a 
good measure of the real worth of the control systems. 

3. An investigation could be made into the effects of allowing 
rod travel in excess of the steady state limit by some percentage 
in an effort to speed up the system response. 

4. A study of the two loop response could be undertaken including 
the effects of coolant loop isolation or steam stop valve closure. 

5. This study used equal transport delays between the reactor 
and the heat exchanger and vice versa. Other plants have substan- 
tially more delay in the cold leg than the hot leg and simulation 


of this could generate new controls. 





APPENDIX 


THE COMPUTiR SIMULATION 


Initially the reactor plant was simulated on the IBM-360 using 
the INTEG program. | The limit of 4500 points required a delta time 
of greater than 0.5 seconds to complete the up power transient. 

In fact, even using a delta time of 0.01 seconds caused over- and 
underflow conditions. The best delta time found was 0.005 seconds 
which only allowed 22.5 seconds of the simulation. 

Discarding the INTEG approach, the IBM Digital Simulation 
Language (DSL) was selected. For all the analysis presented a 
delta time of 0.005 seconds was used in the fixed step size, fourth 
order Runge-Kutta integration. As a check a delta time of 0.001 
seconds was also used which yielded less than 0.3 F difference in 
temperature values and less than 0.4% difference in power levels. 
But since the smaller delta time required over forty-five minutes 
of computer time vice the nine minutes for the larger delta time, 
the 0.005 delta time was selected to increase turn around time. 

The fixed step size integration method was chosen to facilitate 
the calculation of control rod worth. A NOSORT/SORT block was 
required for the control system which precluded the use of an 
integration statement for CRW. Thus a fixed step size was needed 
to be able to generate a constant withdrawal rate. when initially 
running the simulation the actual rate observed was five times as 


‘Inerc is a locally developed integration and plot pvackage of 
the w. R. Church Computer Center, Naval Postgraduate School, Mon- 
terey, California. 
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fast as calculated. As no explanation could be found, other than 
@ possible peculiarity of the program method, rod speeds of two 
and four inches per minute were put into the program resulting in 
ten and twenty inches per minute for the simulation. 

An inconvenience with DSL was that the graphical output on the 
CALCOMP plotter at the Naval Postgraduate School placed the labeling 
for the abscissa to the right of the axis into the graph proper. 
This required modification of the plots for presentation as figures 


in this paper. 
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COMPUTER PROGRAM 


DOWN POWER EXAMPLE 


TITLE DOWN POWER TRANSIENT W/DISCRETE FLOW CHANGE 
CONTRL DELT=0.005 , FINTIM=450 
INTEG RKSFX 
INCON TSIC=458 ,TCCIC=502 ,TCAVIC=508 ,TFIC=1108,THIC=514 
INCON CRW=0,CRP=0,NIC#=1,N=1,CIC=1,C=1,M=2, DLTP=0.05 
NOSORT 
* POWER TRANSIENT INPUT 

EP uUVeS ies) GO Oa 

ie iMac h 410) <eOatOu2 

PSI =1-(TIME-50)*0.08/6 


GO TO 4 
1 PST=1 

COPLee> 
2 PSI=0.2 
3 CONTINUE 


* FLOW CHANGE CONTROL 
cite Oss CO. LO ay 
M=1 

7 CONTINUE 

SORT 
P=0.02*(TAV-TS ) 
TSDOT=70* (P-PSI) /3/M 
TS=INTGRL(TSIC ,TSDOT ) 
TC=((2-0.48/M)*TH+0. 96*TS/M)/(2+0.48/M) 
TCCDOT=(TC-TCC )*M/6 
Tee (NCH ETC Chey Ce Den.) 
TCADOT=(TF-~TCaV-50*M*(TCAV-TCC ) )/25 
TCAV=INTGRL(TCAVIC ,TCADOT ) 
TFDOT=(600*N-TF+TCAV )/25 
TF=INTGRL(TFIC ,TFDOT ) 
THDOT =(2*TCAV-TCC-TH )*M/6 
TH=INTGRL(THIC , THDOT ) 
TAV=(TH+TC )/2 
RLMT=0.0138* ( PSI-1 ) 

NOSORT 

* ROD MOTION CONTROL 
IF ((508-TAV).GT.3) GO TO 4 
Te (0500-PAV oi. = 5) GOTO 5 
IF C(W-PSL).GT.DLTP) GO TO 5 


GO 1026 
4 CRW=CRW+78-7/3 
GOnrO.G 
5 ie (CRW. LE REE) CO.10 6 


Reel. 0.15) GO TO-o 
CRW=CRW-7-7/6 

6 CONTINUE 

SORT 


1) 





KEP=1.120988-1.88E-4*TCAV-2. 3E-5* TF+CRW 
NDOT =2E4* (N* (KiP-1 ) /KEF+7558-4* (C-N) ) 
N=INTGRL(NIC ,NDOT) 
CDOT=0.1*(N-C ) 
CatNT Chin CiG,.c nol, 
CRP=CRW/7E-4 

BeInr 2,97H,TAV,TC.To,1.rol,e,GRP 

END 

STOP 
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